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Abstract 

The magnetic flux penetration dynamics of type-II superconductors in the flux flow regime is studied 
by analytically solving the nonlinear diffusion equation for the magnetic flux induction, assuming 
that an applied field parallel to the surface of the sample and using a power-law dependence of 
the differential resistivity on the magnetic field induction. An exact solution of nonlinear diffusion 
equation for the magnetic induction is obtained using a well known self-similar technique. 
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§1. Introduction 

Theoretical investigations of the magnetic flux penetration dynam- 
ics into superconductors in a various regimes with a various current- 
voltage characteristics is one of key problems of electrodynamics of 
superconductors. Mathematical problem of theoretical study the dy- 
namics of evolution and penetration of magnetic flux into the sample 
in the viscous flux flow regime can be formulated on the basis of a 
nonlinear diffusion-like equation [1-3] for the magnetic field induction 
in a superconductor [4-14], The dynamics of space-time evolution of 
the magnetic flux penetration into type-II superconductors, where the 
flux lines are parallel to the surface of the sample for the viscous flux 
flow regime with a nonlinear relationship between the field and current 
density in type II superconductors has been studied by many authors 
[4-7]. The magnetic flux penetration problem was theoretically stud- 
ied for the particular case, where the flux flow resistivity independent 
of the magnetic field by authors [5]. Similar problem has been con- 
sidered in [6] for the semi-infinite sample in parallel geometry. The 
situation, where flux flow resistivity depends linearly on the magnetic 
field induction was considered analytically in [4]. Analogical problem 
for the creep regime with a nonlinear relationship between the current 
and field has been considered in [8-14], The magnetic flux penetration 
into the superconductor sample, where the flux lines are perpendicular 
to the surface of the sample is described by a non-local nonlinear dif- 
fusion equation [7]. This problem has been exactly solved by Briksin 
and Dorogovstev [7] for the case thin film geometry in the flux flow 
regime of a type-II superconductors. 

§2. Objectives 

In the present paper the magnetic flux penetration dynamics of 
type-II superconductors in the flux flow regime is studied by analyt- 
ically solving the nonlinear diffusion equation for the magnetic flux 
induction, assuming that an applied field parallel to the surface of the 
sample and using a power-law dependence of the differential resistivity 
on the magnetic field induction. An exact solution of nonlinear diffu- 
sion equation for the magnetic induction B(r, t) is obtained by using a 
well known self-similar technique. We study the problem in the frame- 
work of a macroscopic approach, in which all lengths scales are larger 
than the flux-line spacing; thus, the superconductor is considered as 
an uniform medium. 



§3. Formulation of the problem 

Bean [15] has proposed the critical state model which is success- 
fully used to describe magnetic properties of type II superconductors. 
According to this model, the distribution of the magnetic flux density 
B and the transport current density j inside a superconductor is given 
by a solution of the equation 

rotB = p J. (1) 

When the penetrated magnetic flux changes with time, an electric field 
E(r, t) is generated inside the sample according to Faraday's law 



rotE -- 



dB_ 
dt ' 



(2) 



In the flux flow regime the electric field E(r, t) induced by the moving 
vortices is related with the local current density j(r, t) by the nonlinear 
Ohm's law 

E = pj. (3) 
In combining the relation (3) with Maxwell's equation (2), we obtain 
a nonlinear diffusion equation for the magnetic flux induction B(r,t) 
in the following form 



— V \p(B)VB 



(4) 



dB 

dt hq 

Formally, this differential equation is simply a nonlinear diffusion equa- 
tion with a diffusion coefficient depending on magnetic induction B. 
The parabolic type diffusion equation (4) allows to obtain a time and 
space distribution of the magnetic induction profile in a superconduc- 
tor sample. It is evident that the space-time structure of the solution 
of the diffusion equation (4) is determined by the character of depen- 
dence of the differential resistivity coefficient p on the magnetic field 
induction B. Usually, in real experimental situation [16], the differen- 
tial resistivity p grows with an increase of magnetic field induction 



B4>0 
„2 



— , (5) 

where p n is the differential resistivity in the normal state; i) is the 
viscous coefficient, <fro = irhc/2e is the magnetic flux quantum, H c 2 
is the upper critical field of superconductor [16]. In the case, when 
the differential resistivity p is a linear function of the magnetic field 
induction B an exact solution of the diffusion equation (4) can be 
easily obtained by using the well-known scaling methods [1, 2]. For 
the complex dependence of p(B) it can be use by empirical power-law 
dependence p(B) = B n , where n is the positive constant parameter. 



§4. Basic equations 

We formulate the general equation governing the dynamics of the 
magnetic field induction in a superconductor sample. We study the 
evolution of the magnetic penetration process in a simple geometry - 
superconducting semi-infinitive sample x > 0. We assume that the ex- 
ternal magnetic field induction B e is parallel to the z-axis. When the 
magnetic field with the flux density B e (t) is applied in the direction 
of the z-axis, the transport current j(r,t) and the electric field E(r,i) 
are induced inside the slab along the y-axis. For this geometry, the 
spatial and temporal evolution of magnetic field induction B(r, t) is de- 
scribed by the following nonlinear diffusion equation in the generalized 
dimensionless form [7] 



db 


d 




'db' 








b n 






dt ~ 


dx 




dx _ 





(6) 
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where we have introduced the dimensionless parameters b = B/B e , 
j = j/jc, x' = x/xq, t' = t/r and variables xq = B e /fiQj c is the mag- 
netic field penetration depth in a Bean model; r = p n j 2 fio/B 2 is the 
relaxation diffusion time; q is the positive constant parameter. 

The diffusion equation (6) can be integrated analytically subject 
to appropriate initial and boundary conditions in the center of the 
sample and on the sample's edges. We consider the case, when the 
magnetic field applied to sample increases with time according to a 
power law with the exponent of a > 

b(0,t)=b t a (7) 

Boundary condition (7) is equivalent to a linear increase in the mag- 
netic field with time, which corresponds to a real experimental situa- 
tion. As can be easily seen that the case a = describes a constant 
applied magnetic field at the surface of the sample, while the case 
a = 1 corresponds to linearly increasing applied field, respectively. 
The other boundary condition follows from the continuity of the flux 
at the free boundary x = x p 

b(x p ,t) = 0, (8) 

where x p is the dimensionless position of the front of the magnetic 
field. The flux conservation condition for the magnetic field induction 
can be formulated in the following integral form 



/ 



b(x,0)dx = 1. 



(9) 



It should be noted that the nonlinear diffusion equation (6), completed 
by the boundary conditions for magnetic induction, totally determines 
the problem of the space-time distribution of the magnetic flux pen- 
etration into superconductor sample in the flux flow regime with a 
power-law dependence of differential resistivity on the magnetic field 
induction. Solution of this equation gives a complete description of 
the time and space evolution of the magnetic flux in a sample. 

§5. Scaling solution 

In the following analysis we derive an evolution equation for the 
magnetic induction profile and formulate a similarity solution for the 
b(x, t). As can be shown that the nonlinear diffusion equation (6) can 
be solved exactly using well known scaling methods [1, 2]. At long 
times we present a solution of the nonlinear diffusion equation for the 
magnetic induction (6) in the following scaling form 



6(*,*)=t°/(^) 



(10) 



The similarity exponents a and /3 are of primary physical importance 
since the parameter a represents the rate of decay of the magnetic 
induction b(x, t), while the parameter /3 is the rate of spread of the 
space distribution as time goes on. Inserting this scaling form into 
differential equation (6) and comparing powers of t in all terms, we 
get the following relationship for the exponents a and f3 

a + 1 = a(n + q) + ,9(1 + q). 

Using the condition of the flux conservation (9) we obtain 

1 

n + 2q 

which suggests the existence of self-similar solutions in the form 



a = P ■ 



b(z) = t -V(«+2«) z = x ./ t l/(«+2<7). 



(11) 



Substituting this scaling solution (11) into the governing equation (6) 
yields an ordinary differential equation for the scaling function f(z) in 
the form 



(n + 2q) 



dz 



r 



+ 2— +/ : 

dz 



0. 



(12) 



The boundary conditions for the function f(z) now become 

/(0) = 1, /(zo) = 0. (13) 



The above equation (12), depending on the initial and the boundary 
conditions describes a scaling — like behavior magnetic flux front with 
a time — dependent velocity in the sample. After a further integra- 
tion and applying the boundary conditions (13) we get the following 
solution of the problem 



/(*o) 



/(*) = fM 

n + q — 1 



i-(AV +9 



l/(n+9-l) 



(14) 



1 + q \n + 2q 



1 



!/«' 



q/(n + q-l) 



Z Q 



The position of the front zo can now be found by substituting the 
solution (14) into the integral condition (9) and it is given by 



(n+2q)/(n + q-l) 
2 



1 



n + 2q 



1/9 



q/(n+q-l) 



(15) 



n + q 



9+1 



It is convenient to write the self-similar solution (14) in terms of a 
primitive variables, as 



b(x,t)=bo(t) 



(16) 



b (t) = t-V(«+2?) 



+ 9-1 
1 + 9 



^0 

n + 2q 



q/(n+q-l) 



This solution describes the propagation of the magnetic field into the 
sample, the magnetic induction being localized in the domain between 
the surface x=0 and the flux front x p . This solution is positive in the 
plane x p > x 2 and is zero outside of it. Note, that only the x > 
and t > quarter of the plane is presented, because of it has physical 
relevance. The penetrating flux front position x = x p (t) as a function 
of time can be described by the relation 



dx p 
dt 



t -(2«+n-l)/(n+2 9 ) _ 




Curve 1: n=3 
Curve 2: n=7 
Curve 3: n=11 



Fig.l. The profile of the magnetic flux front velocity at different 
values of n=3, 7, 11. 

The velocity of the magnetic flux front decreases rapidly as the 
magnetic flux propagates (Figl). 
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§6. Particular case 

The spatial and temporal profiles of magnetic flux penetration in 
the sample depends on the set of three independent parameters, n, q 
and a. It is of interest to consider the nonlinear diffusion equation for 
the magnetic induction for different values of the exponents n, q and 
a. For a given parameter set n, q and a the form of the scaling func- 
tion f(z) can be obtained by solving the nonlinear diffusion equation 
(6) analytically by a self-similar technique. We solve the nonlinear dif- 
fusion equation analytically to provide expressions for the time-space 
evolution of the magnetic induction for different values of exponents n, 
q and a. Next, we systematically analyze the effect of different values 
of exponents on the shape of the magnetic flux front in the sample. 
Varying the parameters of the equation, we may observe a various 
shapes of the magnetic flux front in the sample. A similar approach 
has been presented in Ref. [7] within the framework of non-linear flux 
diffusion in transverse geometry. As can be shown below that different 
values exponents n and q generate different space-time magnetic flux 
fronts. Below we consider a few more practically relevant examples for 
which the magnetic flux front has a different shape depending on the 
different values of exponents n and q. 

§6. 1. Case q = 1 

Let us first consider the most interesting case q=l. In this par- 
ticular case the spatial and temporal evolution of the magnetic flux 
induction is totally determined by the parameters n and a. In the 
following analysis we derive an evolution equation for the magnetic 
induction profile and apply the scalings of the previous section to for- 
mulate a similarity solution for the b(x, t). For this particular case 
nonlinear diffusion equation (6) can be solved exactly using the scaling 
method. Thus, based on the scalings described in the previous section, 
we get the following relation for the exponents 



1 



The last relation suggests the existence of solution to equation (6) of 
the form 

b(x,t) = t- 1 / ( - n+2 )f(z), z = x/t 1 ^ n+2 >. (17) 
Substituting the similarity solution (17) into the governing equation 
(6) yields an ordinary differential equation for the scaling function f(z) 



dz V rfz / dz 



0. 



(18) 



Integrating the equation (18) by parts and applying the boundary con- 
ditions (13) give 



/(*) 



l/r, 



l/n 



(19) 



_2(n + 2) 

which is the explicit form of the similarity solution, which we have been 
seeking. The position of the front zq can now be found by substituting 
the last solution into the integral condition (9), so we have 



l/n 







1 



_2(n + 2) 

By using the following transformation 



l/n 



dz 



20 sincj, 



and after integrating we obtain 



fn+2)/n 
2 



2(n + 2) 



l/r. 



r i 



It is convenient to write the self-similar solution (19) in terms of a 
primitive variables, as 

,-| l/n 



b(x, t) = 6q 



1 - 



(20) 



where 



bo 



-i l/n 



t -l/(n+2)_ 



.2(n + 2) 

Equation (20) constitutes an exact solution of the nonlinear flux- 
diffusion equation for the situation, when q=l. As can be seen the 
solution (20) describes the propagation of the flux profile inside the 
sample. The profile of the the normalized flux density b(x, t) for this 
case is shown schematically in figure 2. 



b(l.t) 




Fig. 2. The distribution of the normalized flux density b(x,t) at 
different times t=0.1, 0.2, 0.3 for n=l, q=l. 

The penetrating flux front position x = x p (t) as a function of time 
can be described by the relation 

Xp = Zot V(n+2), 

The velocity of penetration of a magnetic flux into a superconductor 
can be naturally determined from the last relation 

V ~ t -<,n+X)/(n+2) 

Interestingly, that the normalized current density j(x,t) in the region, 
< x < x p can be presented using the equations 



1 db 

c dx 



After a simple analytical calculation, we can easily obtain the space 
and time profiles of the normalized current density j(x,t) in the fol- 
lowing form 



26o 



1/n-l 



§6. 2. Case n = 

Let us consider the case n = 0. For this particular case the long- 
time asymptotic behavior of the magnetic induction has the scaling 
form 



b(x,t)=f(z), Z: 



■ xt 



-1/(9 + 1) 



(21) 



Substituting the scaling solution into the governing equation (6) 
yields an ordinary differential equation for for the distribution f(z) in 
the form 



(9 + 1) 



dz 



q + z df 
dz 



0. 



Integrating the last equation and applying the boundary conditions 
(13) we get the following solution of the problem 



3 



3{z) = - 



dz 



-S4-* 2 ) 



1/(9-1) 



(22) 



[29(9 + 1) 

where the position of the front zrj can be found by substituting the 
solution (22) into the integral condition (9), and it is given by the 
following expression 



zo ■ 



2q 



9+1 
<Z-1 



9-1 



2 r( 3 - + - 

\2 q 



1-9' 



1/(9+1) 



The scaling solution (22) can be written in terms of the primitive vari- 
ables as 



1 



£1/(9+1) 



9-1 



1 V(9-l) 



L2 9 (9 + l) 



,2/(9-1) 
'0 



1/(9-1) 



The the flux front can be approximately given as x p = zot 1 /( q+1 \ 
The velocity of penetration of a current density into a superconductor 
is determined from the relation 

v = d ^- ~ t -9/(9+l) 

dt 

We note that, an analogous problem has also been studied in [10, 11] 
in connection with the magnetic relaxation of a superconducting slab 
in the flux creep regime in the framework of an approximate power-law 
dependence of the electric field E on the current density j. The authors 
[11] showed that in the case logarithmic barriers the relaxation process 
causes the system self-organize into critical state. Supposing that the 
homogeneous magnetic induction Bo is induced by a constant mag- 
netic field, they found the expression for the magnetization moment 
in the limit of n > 1. It has been shown by Koziol [12] that a pinning 
potential depending logarithmicaly on current density leads to a sim- 
ilar nonlinear diffusion equation for the spatiotemporal evolution of 
the flux density with a a power-law current-voltage characteristic. An 
approximate and exact solutions of the diffusion problem have been 
derived assuming that an external magnetic field directed to parallel 
of the surface of a sample. The scaling relation between the character- 
istic relaxation time, magnetic field and a sample size has been found. 
Similar problem has been considered by Gilchrist [8, 9]. The flux creep 
problem has been solved by Wang et.all., [13] and for an exponential 
model by authors [14], numerically. 



§6. 3. Case n = 1 

Let us now consider the case n=l. In this particular case the 
spatial and temporal evolution of the magnetic flux induction is de- 
termined by the parameters q and a. In the following analysis we 
derive an evolution equation for the magnetic induction profile for the 
case n=l and apply the scalings of the previous section to formulate 
a similarity solution for the b(x, t). For this particular case nonlinear 
diffusion equation (6) can be solved exactly using the scaling method. 
Thus, based on the scalings described in the previous section, we get 
the following relation for the exponents 

1 

a = (3 



29 + 1 

The last relation suggests the existence of solutions of the form 

b(x,t) = t- 1 /^+ 1 ^f(z), z = x/t 1 ^ 2 "+ 1 l (23) 

Substituting the similarity solution (23) into the governing equa- 
tion (6) yields an ordinary differential equation for the function f(z) 



dz 



(24) 



Integrating the equation (24) by parts and applying the boundary con- 
ditions give 



/(*) 



29 + 1 



1/9 



1 - 



20 



(9 + l)/9 



(25) 



which is the explicit form of the similarity solution we have been seek- 
ing. The position of the front can now be found by substituting the 
last solution into the integral condition (9), and we have 

Z0 = [ g l/(9+l)(2g+l)])('?+ 1 )/(^+ 1 ). 

It is convenient to write the self-similar solution (25) in terms of a 
primitive variables, as 



where 



b(x,t) 



b(x,t) 



-1/(29 + 1) 



1/9 



(26) 



29 + 1 



1/9 



Equation (26) describes an exact solution of the nonlinear flux- 
diffusion equation for the situation, when n=l. As can be seen the 
solution (26) describes the propagation of the flux profile into the sam- 
ple. The the flux front can be approximately given as x p = ti/( 2 9+i) 
The velocity of penetration of a magnetic flux induction front into a 
superconductor is determined from relation 

v= dXp ^ t -2,/(2,+ l) 
dt 

The velocity of the magnetic flux front decreases rapidly as the mag- 
netic flux propagates. 



§6. 4. Case n = 1, 9 = 1 

Let us now consider the case, where n = 1, 9 = 1. For this 
particular case the nonlinear diffusion equation for the distribution 
magnetic flux admits an exact self-similar solution with a sharp front 
moving with a constant velocity. The scaling analysis show that the 
diffusion equation can be solved for the following values of parameters 
1). a=0, 13=1/2; 2). a=l/3, (3=2/3; 3). a=l, /3=1; 4). a=l/3, 
(3=1/3. Let us consider a solution of the nonlinear diffusion problem 
for the exponents a=l/3, (3=1/3. In this case the nonlinear diffusion 
equation (6) admits an exact solution [4] with the similarity variable 

b(x,t) = t- 1 / 3 f(xt- 1 ^). (27) 

By substituting the solution (27) into diffusion equation (6) we obtain 
an ordinary differential equation for the scaling function f(z) in the 
form 



dz 



('£) 



dz 



(28) 



Integrating the equation (28) and using the boundary conditions for 
the magnetic induction f(z) we get a simple equation 



3— + z = 0. 

dz 



(29) 



Separating variables and integrating equation (29) gives an implicit 
solution in the form 



(30) 



where the integration constant is determined by integral relation (9) 
and has the form 



20 



1/3 



Thus, we have solution for the magnetic induction in a primitive vari- 
ables 



4 



t) 



6tV3 



9 \2/3 
2, 



x 2 
t 2 / 3 



(31) 



The expression (31) describes an exact solution of the nonlinear flux- 
diffusion equation for the situation, when n = 1, q = 1. As can be 
seen the solution (31) describes the propagation of the flux profile into 
the sample. The profiles of the the normalized flux density b(x, t) and 
electric field density e(x, t) for this case is shown schematically in the 
figures 3a-3b. 



b(K,t) 




e(*,t) 



Curve 1. t=0.1 
Curve 2: t=0.2 
Curve 3: t=D.3 




Curve 1: t=0.1 
Curve 2: t=0 2 
Curve 3: t=0.3 



Finally, in the case q=0, /3=l/2 the flux front position has the 
form of x p oc t 1 / 2 and an approximate analytical solution for the case 
can be easily found solving the diffusion equation with the help of sim- 
ilarity variable b(x,t) = f(z). where z = /(x/i 1 / 2 ). In the vicinity of 
the front (x oc x p ) the magnetic flux profile is given by the following 
simple relation b oc \/2{x — x p ). In this case, which corresponds to a 
constant external magnetic field the flux front position x p grows with 
time as x v oc 4 1 / 2 for long times. 



§6. 5. Case n = 2, 



The scaling analysis show that the diffusion equation can be solved 
for the following values of parameters l)a=l/4, /3=l/4; 2) a=l/4, 
/3=l/6. The solution to diffusion equation for this case can be pre- 
sented in the following form 



b(x,t) = t- 1 ^ 4 f(x/t 1 ^ 4 ). 



(34) 



Substituting the solution (34) into the partial differential equation (6) 
and integrating with the boundary conditions we can see that magnetic 
flux profile has an exact solution in the form 



' 2 



z 2 

X ~V 2 

z o 



1/12 



where the integration constant has the form 

2 

20 = —p=. 



The solution for the magnetic induction in a primitive variables has 
the form 



b(x, t) 



1 



4 tV2 



1/2 



(35) 



The expression (35) describes an exact solution of the nonlinear flux- 
diffusion equation for the situation, when n = 2, q = 1. As can be 
seen the solution (35) describes the propagation of the flux profile into 
the sample. The evolution of the self-simulating process of magnetic 
field penetration into a superconductor is shown schematically in the 
figures 4a-4b. 



Figs. 3a and 3b. The distributions of the normalized flux density 
b(x,t) and electric field density e(x,t) at different times t=0.1, 0.2, 
0.3 for n=l, q=l. 

It is evident from the solution (31) that the shock wave like mag- 
netic flux propagates with wave front 

x p = z * 1/3 . (32) 
It can be shown that the sharp flux front moves with the time 
dependent velocity 



dx p 
dt 



z t 



-2/3 



(33) 



The case q=1/3, /3=l/3 implies that the flux front asymptotically ex- 
pands according to power-law x p oc t 1 / 3 and magnetic flux penetration 
profile decreases as b oc t — 1 / 3 . Similar problem for the case penetra- 
tion of magnetic field into superconductors has been studied by Bass 
[4]. An analogous theoretical result in the viscous flow mode for vor- 
tices was obtained in [5] using the model of the critical state for oxide 
high-T c superconductors. 




Curve 1: t=0 1 
Curve 2- t=0 2 
Curve 3: t=0 3 
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Curve 1: t=D.1 

Curve 2: t=0.2 

' Curve 3: 1=0.3 

Figs. 4a and 4b. The distributions of the normalized flux density 
b(x,t) and electric field density e(x,t) at different times t=0.1, 0.2, 
0.3 for n = 2, q = 1. 

In this case the front coordinate is x p = i 1 / 4 . A main result is that 
the magnetic flux will move as a shock front with velocity v = t 3 / 4 . 
Similar problem for the case penetration of magnetic field into super- 
conductors has been studied in [5] . 



§6. 6. Case n = 0, q = 1 

In this case the diffusion equation can be solved in a different val- 
ues of a and j3, in particular 1). a=0, /3=l/2; 2). o=-l/2, £=-1/2; 
3). a=X, /3=1. Let us consider the case o=0, /3=-l/2. In this case the 
differential equation is linear and admits an exact analytical solution. 
At long times we present a solution of the nonlinear diffusion equation 
for the magnetic induction (6) in the following scaling form 



b(x,t)=t a f[^\. (36) 
Inserting the scaling form (36) into differential equation (6) we have 

e+^ll=r-^Lf + ^fY (37) 

Comparing powers of t in all terms in (37), we get the following rela- 
tionship for the exponents a and /3 



0: 



Then the problem admits a similarity solution of the form 

b = t- 1 / 2 f(xt- 1 / 2 ). (38) 
Substituting the expression (38) into the nonlinear diffusion equation 
(6) gives an ordinary differential equation for the function similarity 
function f(z), namely 



d?f 

dz 2 



+ z— + f - 
dz 



0. 



(39) 



Integrating (39) twice and applying the boundary conditions we have 
a self-similar solution in the form 



b(x,t) 



where 



: exp 



/it. 



(40) 



The last formulae describes an exact solution of the nonlinear flux- 
diffusion equation for the situation, when n = 2, q = 1. As can be 
seen this solution describes the propagation of the flux profile into the 
sample. The profile of the self-simulating process of magnetic field 
penetration into a superconductor for this case is shown schematically 
in the figures 5a-5b. 



Curve 1 t=0.1 
Curve 2: t=0 2 
Curved. 1=0 3 
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Curve 1: 1=0.1 
Curve 2: 1=0.2 
Curve 3: t=0.3 



Figs. 5a and 5b. The distributions of the normalized flux density 
b(x,t) and electric field density e(x,t) at different times t=0.1, 0.2, 
0.3 for n = 0, q = 1 . 

Similar problem has been solved in [6] using an extended critical 
state model for type-II superconductors. Assuming that the critical 
current density and resistivity are constants a one-dimensional diffu- 
sion equation for the magnetic field were solved analytically for the 
case when an external magnetic field is increased with time according 
to a power law with the exponent of one- half. 

§7. Flux creep 

Let us consider the magnetic flux penetration process for the flux creep 
regime. According to Kim-Anderson theory [17, 18] the thermally 
activated flux motion in superconductor sample is described by an 
Arrhenius-type expression 

v = v exp[-U/kT], (41) 

where vo is the resistivity at T=0, U is an activation energy for ther- 
mally activated flux jumps, and determines the vortex pinning; T is the 
temperature and k is the Boltzmann constant. The activation energy 
U = U(j, B,T) depends on temperature T, magnetic field induction 
B and current density j. For the simple case it can be presented by 
Kim-Anderson [17] formulae 



U(j) = Uo[l 



(42) 



here Uo is the characteristic scale of the activation energy and j c = 
jc(B) is the critical current density. In the flux creep state the effective 



G 



activation energy U grows logarithmically [19] with decreasing current 
density as 



U(j) = U 



(43) 



where the exponent n depends upon the flux creep regime. The ex- 
pression (43) gives a quite realistic description for activation barriers 
in a wide range of temperatures and magnetic fields. The power law 
characteristic dependence for U(j) has been observed in numerous ex- 
periments, and it has been used extensively in recent theoretical stud- 
ies for the field and current distributions in superconductors [19, 20]. 
Equation (43) is equivalent to a power law for the flux-flow resistivity 
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(44) 



In this case the phenomenological relation E(j) may be chosen in the 
power-law form 
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(45) 



If n = 1 the last equation reduces to Ohm's law, describing the nor- 
mal or flux-flow state. For infinitely large n, the equation describes 
the Bean critical state model j = j c [15]. When 1 < n < oo, the last 
equation describes nonlinear flux creep. 

In such case the analytic solution of the nonlinear creep equation 
can be constructed by choosing the critical current dependence on 
the magnetic field. Many models have been for the functional form 
of j c (B). For the critical current we adopt the power- law model [21], 
which can be applied over a relatively wide magnetic field range except 
in the high field region near the upper critical field 



jc(B)=j a 



'Bo\ 



(46) 



where jo and Bo are the characteristic values of the current density 
and magnetic field induction; 7 is the dimcnsionlcss pinning parameter, 
usually < 7 < 1. If we assume 7=0, the above model reduces to the 
Bean-London model [15]. This model is applicable to the case where 
j c can be regarded approximately field independent. This power-law 
model for the critical current has also been used by other groups [14]. 
Another possible decay law would be exponential, which has often been 
used to take into account its decrease with the magnetic field [22] 
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where Bo is a phenomenological parameter related to the pinning abil- 
ity: the smaller it is, the more drastic is the decrease of the critical 
current with field. The numerical methods have been applied to re- 
solve the flux diffusion equation, employing the exponential critical 
state model [14]. 

Next, based on the power law and exponential models, we shall 
study the distribution of the magnetic induction, current density and 
magnetization of superconductors. 

§7.1. Power law model 

For the one-dimensional geometry the spatial and temporal evolution 
of magnetic field induction B(r, t) is described by the following non- 
linear diffusion equation [10] in the dimcnsionlcss form 



db _ d 

dt dx 
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n + l 



db 

dx 



(47) 



where we have introduced the dimcnsionlcss variables and parameters 
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The solution of above parabolic type diffusion equation allows to ob- 
tain the time and space distribution of the magnetic induction profile 
in the considered sample. Inserting the scaling form (10) into differen- 
tial equation (47) and comparing powers of t in all terms, we get the 
following relationship for the exponents a and P 

a + 1 = p + a ( 7 rt + 1) + an + 0n. (48) 
Using the condition of the flux conservation we obtain 

1 



(2n + n-y + 1) 

which suggests the existence of self-similar solutions in the form 

b(x,t) = t- 1 /( 2n+n T+ 1 )/(z), z = x/t 1 ' {2n+n ' 1+1 \ (49) 
Substituting this scaling solution (49) into the governing equation (47) 
yields an ordinary differential equation for the scaling function f(z) in 
the form 
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(50) 
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The solution of the equation (50) must satisfy the following boundary 
conditions for the scaling function f(z) 

/(0) = 1, f(zo) = 0. (51) 
After further integrating of equation (50) and applying the boundary 
conditions we get the following solution of the problem 
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(52) 
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The position of the front zo can now be found by substituting the 
solution (52) into the integral condition (9) and it is given by 



7 + 2 1 
7 + 1 2 



n + 



(2n+n 7 +l)/( 7 +l) 
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The last solution can be presented as 
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b(x,t) = b 

where 

6 (0, t) = t -l/(2n+»7+l) 
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(53) 
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This solution describes the propagation of the magnetic field into the 
sample, the magnetic induction being localized in the domain between 
the surface x=0 and the flux front x p . The flux front can be approxi- 
mately given as x p = sot 1 /' 2n + n7 + 1 ' . The velocity of penetration of 
the magnetic flux induction front into the superconductor sample is 
determined from the relation 

V ~ t -«(2 + 7)/(2n+n 7 + l)^ ^ 

The velocity of the magnetic flux front decreases rapidly as the mag- 
netic flux propagates. 

Let us consider the most interesting case n=l. In this particular 
case the spatial and temporal evolution of the magnetic flux induction 
is totally determined by the parameters 7, a and p. In the following 
analysis we derive an evolution equation for the magnetic induction 
profile for the case n=l and apply the scalings of the previous section 
to formulate a similarity solution for the b(x, t). Based on the scalings 



7 



described in the previous section, we get the following relation for the 
exponents 



a = /3 



1 



7 + 3 

The last relation suggests the existence of solutions of the form 
b(x,t) = t-VC-r+a)/^), z = a! / t i/(-r+3)_ 



(55) 
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Substituting the similarity function (56) into the governing equation 
(47) yields the following solution of the problem 
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which is the explicit form of the similarity solution we have been seek- 
ing. The position of the front has the form 
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Fig. 6a and 6b. The distributions of the normalized flux density 
b(x,i) at time t=l, for 7=1, 2. 
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(58) 



(59) 



where 
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t -l/(7+3) 



Equation (59) constitutes an exact solution of the nonlinear flux- 
diffusion equation for the situation, when n=l. The evolution of the 
self-simulating process of magnetic field penetration into a supercon- 
ductor is shown schematically in figure 6a and 6b. 




§7.2. Exponential model n=l 

Let us consider a solution to the nonlinear diffusion equation for 
the exponential model, assuming that n=l. The nonlinear differential 
equation (47), describing evolution of the magnetic flux penetrated 
into the superconductor sample for the exponential model can be trans- 
formed to the following form 



, db 



dx 



(60) 



db _ d 
dt dx 

It is remarkable that one of explicit solution of the nonlinear diffusion 
equation (60) can be obtained by using the method of separation of 
variables. According to the method of separation of variables, we look 
for the solution of (60) in the form 

b(x,t) = a(x)g(t). (61) 
Substituting this expression for b(x, t) into the partial differential 
equation (60) and separating variables we obtain the following equa- 
tions for the distributions of the functions g(t) and a(x) 



dt 
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Integrating the differential equations (62) and (63) with respect to t 
and x, respectively and using the boundary conditions we obtain 
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(64) 
(65) 
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where t p is a constant, the parameter x p is defined by using the inte- 
gral relation (9). Combining the equations (64) and (65) according to 
relation (61) we get 
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§8. Magnetization 

Let us consider the magnetization profile of the sample which is 
given by the following relation 



-HaM{t) = b 
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Substituting here the scaling law, Eq. (15) one obtains 

Xp(t) 



-HoM(t) = b 
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where d is the distance of the flux front penetrating into the sample; 
4> n is a numerical factor which is determined by the average value of 
the scaling function f(z) in the region < x p < d, so 
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Expressing the time dependence explicitly, the magnetization can be 
written as 



-HoM(t) = b 



1 — Ipn 



for < t < t* , where t* = {d/zo)^ ■ Differentiating the last equation 
and using an explicit solution for x p (t) for the case 7 = — 1/n one 
obtains the expression to the magnetic relaxation rate 



dM 
d\nt 
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n + 1 



1/(1+1) 



As the parameter n increases the profile of the relaxation rate is seen 
to become more linear. At low temperature limit when n S> 1 the 
relaxation rate varies linearly with time. The time dependence of the 
magnetic relaxation rate is shown in figure 7. 
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Curve 1: n=3 
Curve 2: n=5 
Curve 3; n=11 



Fig. 7. The time dependence of the magnetic relaxation rate for n=3, 

5, 11. 

It has been shown by Vinokur [11] that at low temperatures the 
relationship between the time and relaxation rate is a linear over a 
very long time period in the early stage of the flux penetration. A 
such perfect linear relationship characterizes that the system behaves 
a self-organized criticality. 



Conclusion 

Thus, the problem of magnetic flux penetration into the half-space 
superconductor sample is studied in the flux flow regime in parallel ge- 
ometry assuming that an external magnetic field increasing with time 
in accordance with the power law (7). Assuming that the flux flow 
resistivity as a power-law function of the magnetic field induction we 
found an exact analytical solution for the nonlinear local magnetic flux 
diffusion equation. The obtained solution describes space-time distri- 
bution of the magnetic induction in the sample. It was shown that 



the magnetic flux density profiles at different times follow the different 
scaling law. Similar scaling profiles can be found for the electric field 
and current density profiles, too. It was also shown that the spatial and 
temporal profiles of magnetic flux penetration in the sample depends 
on the set of three independent parameters, n, q and a. For a given 
parameter set n, q and a the form of the scaling function b(x, t) was 
obtained by solving the nonlinear diffusion equation for the magnetic 
field induction analytically by a self-similar technique. We analyzed 
the effect of different values of exponents on the propagation of the 
flux front in the sample. Varying the parameters of the equation, we 
have observed a various shapes of the magnetic flux front. Finally, we 
have discussed the flux creep problem, briefly. 
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